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The nature of glass phases of hquid crystals in random porous media depends on the effective 
disorder strength. We study how the disorder strength depends on the density of the porous media 
^^ ' and demonstrate that it can increase as the density decreases. We also show that the interaction of 

^.ZJ the liquid crystal with random porous media can destroy long-range order inside the pores. 

o. 

CN ■ PACS numbers: 61.30.Pq, 61.30.Gd 

^^ , In many natural and technological systems liquids are confined in random porous media. Two examples are water 

in soil and oil in porous rocks. Liquid flow in porous media involves complicated dynamical phenomena and this 

\l ' motivates significant research efforts. The physics becomes even richer when a porous matrix confines a liquid crystal. 

Interaction with random media affects not only dynamic but also static properties of liquid crystals. Long-range 

order in liquid crystalline states is unstable in the presence of even arbitrarily weak quenched disorder [1, 2] and novel 

tj^ . glass phases emerge instead [3]. In particular, slow glassy dynamics was reported in [4]. Another disorder effect is 

O ' the suppression of the isotropic-nematic [5] and nematic-smectic [6] phase transitions. Properties of the glass states 

^ . depend on the relative strength of disorder and elastic forces. When disorder is strong, liquid crystals are expected 

to have a finite correlation length [3], while in the weakly disordered case, quasi-long-range ordered glass phases can 

emerge [3, 7, 8]. One would naively expect that the effective disorder is weaker in low density random porous media 

such as aerosil gels [9], i.e., when the pore size is large. However, recent experiments found no signs of ordering in a 

i-rt ' low density aerosil system [11] while liquid crystalline order was detected in a high density bulk fractal silica gels [12]. 

(~| The results of light scattering experiments with nematics in dilute aerosils [10] are consistent with a faster decay of 

O ' correlation functions at large distances in lower-density aerosil gels. 

, ^, , In this paper we examine how the effective disorder strength at scales greater than the pore size depends on the 

pore size and show that this dependence has nontrivial character. In particular, the disorder strength can grow as 
Ps) \ the pore size increases. There are two reasons for this. First, the effective disorder strength is the ratio of the typical 
K*' ■ random and elastic energies. Since the latter grows slowly as a function of the pore size, the former can dominate for 
Cn \ low aerogel densities. Second, there are long range correlations in the structure of random porous media and these 
correlations can enhance the effective disorder. We also consider the physics at scales smaller than the pore size. At 
such scales, random media cannot be viewed as uniform and the fractal structure affects the confined liquid crystal. 
^^ , We show that the interaction of the liquid crystal with the random pore surface can destroy long-range order in the 
^+ bulk of the pore and calculate the order parameter correlation function for a simple model of random fractal porous 

f^ , media. 

"^ ■ The organization of this article is as follows. First, we address the effective disorder strength in systems of continuous 

2 ' symmetry in porous media in the presence of random field and random anisotropy disorder, the latter being relevant 

^ . for nematic liquid crystals. We also consider the "directed random field" universality class which can be realized 

' ' in random gels based on amino-acids [13]. Second, we investigate the problem of a liquid crystal interacting with a 

random Sierpinsky gasket surface. This is the simplest model of what might occur inside a pore. We find scaling 

behavior for the correlation functions in the pore. 

^ , Most of our results are valid for a general system with a continuous symmetry order parameter. We focus on 

;V ■ disordered nematics and the random field XY model [14, 15], the latter being the simplest model with a continuous 

. ^H ! symmetry group. These models represent two different types of quenched disorder: random anisotropy and random 

j^ field. In nematics antiparallel orientations of the director correspond to the same physical state and hence to the same 

?H disorder energy, i.e. the disorder is random anisotropy. Random field disorder is possible in ferroelectric liquid crystals 

where the energy is different for antiparallel dipole moments. We will show that the role of geometric correlations is 

stronger in the random anisotropy case. 

Nematics can be described by the Hamiltonian [3]: 

i^n = y"dV{i[Xi(V-n(r))2 + if2(n(r).Vxn(r))2 + if3(n(r)xVxn(r))2]- ^ /i„^(r)n„(r)n^(r)} (1) 

where Ki are the Frank moduli, the unit vector n is the director, and the random tensor hap{v) describes the 
interaction with the porous medium. We assume that the preferred orientation of the nematic molecules near the 
surface of the pores is parallel to the surface. We also assume that the nematic order parameter is relatively large in 
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the pores, so that it is appropriate to use the director description, rather than the more general tensor nematic order 
parameter. 

The random field XY model has the following Hamiltonian 

H^y = Jd\['-^ Y, (dc^Hrff -hi{r) COS ^{r)-h2{r) sin cf>{r)Y (2) 

a—x,y.z 

where (j)(r) is the polar angle of the spin, and hi{r) and /i2(r) are the x and y components of the random field 
respectively, each with zero average. Random field disorder is possible in ferroelectric liquid crystals. Since ferroelectric 
liquid crystals have other degrees of freedom in addition to their dipolc moment, the random field XY model does 
not provide a full realistic description of ferroelectric liquid crystals. However, this simple random field model does 
provide useful insight. The model (2) is also relevant to discotics and smectics confined in stretched aerogels [8] but 
only at length scales greater than the pore size. 

The models Eqs. (1) and (2) can be used for length scales both greater and smaller than the pore size. However, 
the probability distribution of the random field is independent of the coordinates only in the former case. For short 
length scales the structure of the porous matrix must be taken into account. The effective model at length scales 
greater than the pore size must be derived from the short length scale model. Such a derivation is difficult and only 
order-of-magnitude estimates for the model parameters have been obtained by fitting experimental data [16]. We 
focus on the disorder strength dependence on the aerogel density. The disorder strength in the long length scale model 
determines the type of ordering at long length scales. For weak quenched disorder, quasi-long-range order is expected 
for nematics and the random field XY model. The ordering disappears for strong quenched disorder. 

Let a be the microscopic cut-off scale for the short length scale model, ^ the pore size, and df the fractal di- 
mension of the aerogel. An Imry Ma type estimate [1] shows that the typical elastic energy at the scale ^ is 
J(^) ^ J(a)(^/a)''/(^/a)^ = J(a)(^/a)''^^, for a d-dimensional system and the disorder energy h{£_) ^ h{a) ^JJ^JoY^ . 
The effective disorder strength at the scale ^ is the ratio of the typical disorder and clastic energies: 

/i,(0 - h{^)J{a)/J{S,) ^ ^'^fn-d+2^ (3) 

This effective disorder strength hj. enters the random contribution to the Hamiltonian of the long length scale model, 
if one keeps the same elastic constant as in the short length scale model. We see that the disorder strength in the 
long length scale model grows as a function of the pore size for df > 2(i — 4 (destroying long-range order in this case) 
while the aerogel density p ^ ^'^f~'^ decreases for all fractal objects. Thus, in d = 3, the disorder strength grows as 
a function of ^ for all d/ > 2, and we note that this condition is satisfied in recent experiments on highly porous 
gels [10, 12, 17]. Thus, increasing the pore size can lead to an increase in the effective disorder strength with an 
accompanying decrease in the aerogel density. In particular, if the disorder is strong for a smaller pore size it is also 
strong for a greater pore size. A simple way to understand the critical fractal dimension df ~ 2d — 4 = 2 for d = 3 is 
based on the following observation: The effective disorder strength is infinite for a system of disconnected pores and 
the fractal dimension of the walls of disconnected pores is always 2 or greater. 

The above estimate neglects correlations of the random field h at different points. This approximation is valid 
for the random field disorder. For the nematic problem Eq. (1), geometric correlations in random media result in 
correlations of the random anisotropy since nematic molecules tend to be parallel to the pore surface. This affects the 
effective disorder at the scale ^. The existence of this effect becomes obvious when the porous matrix is not random 
at length scales shorter than ^. Then the random anisotropy is the same at all points of the random medium inside 
the pore. Hence, h{S,) ^ (,'^'' and hr ^ C'^''^^, where df is the fractal dimension of the porous medium. The effective 
disorder strength then grows as a function of ^ for all df > 1. A similar effect is possible for a random fractal. Indeed, 
as df ^ 1 we expect that random and non-random fractals should exhibit similar properties for the case of random 
anisotropy disorder. For example, consider a random fractal whose construction is based on fractional Brownian 
motion [18]. The coordinates of the A^th segment of the random fractal are given by the equation [18] 

N 

B"{N) = J2^^(n){N-n + l)"-^/^, (4) 

n=l 

where the Hurst exponent H is related to the fractal dimension by df — 1/H [18], and AB(n)is a random Gaussian 
vector with n an integer parameter. A Hurst exponent of 1/2 corresponds to ordinary Brownian motion. We imagine 
that this fractional Brownian walk forms the backbone of a tube-like fractal object, and the nematic molecules are 
anchored along the surface of the tube. To obtain a tractable model, we assume that the surface area of the nth 



segment of the tube is proportional to (AB^(n))^, where AB^(n) — B^(n) — B^(n — 1). The average surface area 
of a segment is finite and hence the structure is weU-defined for H < 1. The total strength of the random anisotropy 
associated with each segment hap{n) ^ AB^ ABp is proportional to its area. An estimate of the total disorder 
strength in a pore of linear dimension ^, enclosing a fractal consisting of rimax segments with the z axis chosen to lie 
along the average director is then given by: 
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where, 



7 = 



_ / 1/2, 



2iJ - 1 = 2/df - 1, 



i?<3/4 
H> 3/4 



(6) 



and the bar denotes an average over all realizations of the random fractal. Thus, as df -^ 1 (i.e., H < 3/4), the 
renormalized disorder strength scales with ^ in a manner different than Eq. (3): 



hriO-e- 



(7) 



One can see that the effective disorder decreases most rapidly as the function of the pore size for df = 4/3. Even in 
this case hr ~ ^~^/^ and hence to decrease the disorder strength by one order of magnitude one needs to increase the 
pore size by several orders of magnitude. At df — I, Eq. (7) predicts no dependence of the renormalized disorder 
strength on the pore size. In the presence of any preferred macroscopic orientation for the segments of the random 
surface of the pore, the effective random anisotropy enhancement is stronger than Eq. (7). We do not expect the 
geometric correlations to play an important role in systems with random fields which can be randomly oriented along 
one of the two directions parallel to the segment axis. 

Geometric correlations are, however, important for the "directed random field" universality class. Let a random 
matrix consist of random polymer chains made of identical monomers. We assume that the monomers are asymmetric 
with non-equivalent ends A and B. The A-end always links the i3-end of the next monomer in the chain. The pore 
diameter equals the chain size. Consider a ferroelectric smectic C* confined in such a random medium. Its dipole 
moment couples linearly to the local random field which is parallel to the network monomers [19] and is directed from 
their A-ends to their B-ends. Smectics C* are chiral and form helical structures [20]. However, the typical period 
of the helical structure is large and at smaller scales the liquid crystal can be considered as a uniform ferroelectric. 
Then the standard Imry-Ma argument can be used. 

For any fractal dimension df, the vector sum of all random fields associated with the polymer of length (.^'"^f is 
h(^) ^ ^, where S, is the typical distance between branching points of the network. Hence, the disorder strength in 
the long length scale model with ultraviolet cut-off ^ is scale independent: hr{£,) ~ h{0/^ — const. 

What happens inside a pore? We need to study the limit ^ ^ oo. Let us first consider the exactly solvable Larkin 
model [2] of a system with continuous symmetry: 
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(8) 



where h(r) is a random field which is nonzero on the random fractal only, and its probability distribution is charac- 
terized by h{ri)h{r2) = AS{ri — r2). We want to calculate the correlation function G'(ri,r2) = (0(ri) — (j){r2))'^. Let 
si and S2 denote the distances from the points ri and r2 to the random fractal. We assume that R = |ri — r2| > 
max(si,S2). At low temperatures, the equilibrium configuration of the field (f) can be found by minimizing the 
Hamiltonian (8): 0(r) == / (i^r'/i(r')/47r|r - r'|. Hence, 
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where the integration extends over the random fractal. The result depends on the fractal dimension: for d/ < 2, 
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for df = 2, G ^ logR/si + logR/s2; for d/ > 2, G ~ R''-f ^. As R ^ 00, the correlation 



function diverges for df > 2. In terms of the "physical" order parameter n — (cos </>, sin (/>) , the correlation function 
limfl^oo(n(ri)n(r2)) = limfl;^oo(cos((/)(ri) — (j){r2))) = lim^^oo exp(— G/2) — > 0. Hence, for such dj, long range order 
is destroyed inside the pore. 

We now consider the more realistic and complicated nonlinear model (2). The random field XY model cannot be 
solved exactly and we will use a renormalization group (RG) procedure. We will develop a double e-cxpansion in 
e = 4 — rf and e/ = 4 — 2rf + df, where d is the dimensions of space and dj the random fractal dimensions. This 
choice of small parameters is motivated by the fact that d = A is the critical dimensions for the bulk random field 
XY model [14, 15], and the Imry-Ma argument of Eq. (3) shows that the random fractal destroys long-range order 
for df > 2d — 4 (cf. the Larkin model, where d/ = 2 is the critical fractal dimensions for d — 3). In contrast to the 
Larkin model, the results depend on the structure of the random fractal. We will limit our discussion to a simple 
model based on the Sierpinsky gasket [21]. We divide the system of linear size Lq into blocks of size Li = BLq , 
with 6 < 1. We assume that disorder is present only in a fraction p = {L^/ LiY^^'^ of the blocks. Then we divide 
each block of size Li into blocks of size L2 = bLi. For each block of size Li that contains segments of the random 
fractal we randomly choose a fraction p of subblocks of size L2, such that only these subblocks contain segments of 
the random fractal. We iterate this procedure and generate a random fractal of dimension df. We implement the 
Wilson shell renormalization group (RG) by changing the cut-off scale from _L„ to Ln_i at each step, to generate the 
one-loop RG equations. For the bulk random field XY model this problem was solved in Ref. [22]. In the case when 
disorder is present on a two-dimensional surface only, the problem was studied in Ref. [23] . We will thus only bricfiy 
discuss the RG procedure which is similar to Refs. [3, 15, 22, 23]. 

Due to the periodicity of the Hamiltonian there is no renormalization of the order parameter. There also is no 
renormalization [22] of the elastic term in (2). Hence, the scaling dimension of the temperature [22] is 2 — e. As for 

disorder, we know that all random contributions of the form h\ cos kcj) + h2 sin kef) are relevant operators in the 
vicinity of the zero-temperature fixed point [14, 15]. After replica averaging this corresponds to the following structure 
of the relevant terms in the Hamiltonian: 

-'a ah 

where a, 6 are replica indices, F is the function describing disorder at the ultra-violet cut-off scale i„, P(r) = 1 
in the subblocks of size L„ containing elements of the random fractal, otherwise P{y) ~ 0. Eq. (10) is ob- 
tained by averaging with respect to disorder inside the subblocks of size L„. We can average over the disor- 
der distribution inside the blocks of size i„_i instead. This will include averaging with respect to the distri- 
bution of P(r) which is nonzero with probability p in different points of the blocks of size L„-i which con- 
tain segments of our fractal. Instead of cxp(/ d'^r^^^ P(r)i^((/)a(r) — (/)f,(r))/2T^) in the replica action wc will 

get Hr[(l -p) -f pexp(P'(r)E„fc^(0a(r) - Uv))/2T^)] = exp{J d''rP'ir)[j:^,pF{Mr) - Mr))/2T^ + EabcAP - 
p^)F{(j)a{r) — (j)b{r)) F{(j) c{r) — (j)d{r))/8T'^ + ...]), where -P'(r) = 1 if r belongs to a block of size L„_i containing 
segments of the fractal. We need to keep only the first term X]a6P-^('?^a(r) — 06(r))/2r^ in the integral in order to 
derive the one-loop RG equation. Indeed, we will see that F ^ Cf at the fixed point. We note that (p^ — p) ^ 
as e,e/ -^ 0. Hence, the second term in the square brackets is 0{e^,e^f) and will not contribute to the one-loop RG 
equation. The same is true for the higher-order terms denoted by the dots. Thus, one can directly repeat the standard 
derivation of the RG equations for the bulk random field XY model. The only modification is related to the factor p 
multiplying F. This factor will change the term linear in F in the one-loop RG equation. We obtain: 

dF{cb)/dlnL = e/F(0) + F"((/))V2 - F"(0)F"(O), (11) 

where the primes refer to derivatives of F with respect to its argument and the factor I/Stt^ has been absorbed by F. 
The fixed point solution of the above equation is known [3, 22]. Due to the symmetry of the problem we need a periodic 
solution with period 2tt. At the interval < (/) < 27r it is given by the equation F{(p) — 27r^e/[l/36 — {4>/2tt)'^{1 — 
0/27r)^]/9. This solution corresponds to e/ > 0. Otherwise F ^ 0. The temperature goes to zero at the fixed point. 
We can now calculate the correlation function G — {4>{ri) — 0(r2))^ — 2 J d'^q\<j>'^ {q)\{l — cosq • (ri — r2))/(27r)'^. At 
each step of the RG procedure it is sufficient to use the quadratic in (j) part of the renormalized replica Hamiltonian 
to calculate the contribution corresponding to 1/Ln > q > 1/L„_i. If both points ri and r2 are at a distance 
s < |ri — r2| from the random fractal, we obtain G = (2e/7r^/9) ln(|ri — r2|/s). The correlation function diverges for 
large |ri — r2|, if e/ > 0. Hence, the random fractal destroys long range order in the pore for df > 2d — 4, in agreement 
with our Imry-Ma estimate above. The above result for the correlation function is valid only if the bare disorder at 
the microscopic cut-off scale a is weak since we did a perturbative calculation. 



We have calculated how the effective disorder in models, Eqs. (1) and (2) depends on the pore size. If the disorder 
is weak in the long length scale model with cutoff length ^ then we expect quasi-long range order [3, 14, 15] in both 
models. If the disorder is strong then only short range order is possible [3, 14, 15]. At strong disorder the perturbative 
RG analysis is no longer valid and quasi-long-range order is impossible. The condition of weak disorder means two 
things: (1) the effective disorder at the pore size scale ^ is weak; (2) the renormalizcd disorder is weak at long length 
scales, i.e. e = 4 — d is small. The critical value of the disorder strength is non-universal and depends on microscopic 
details. 

In conclusion, we have found the dependence of the effective disorder in random fractal media on the pore size and 
the fractal dimensions. The dependence is non-trivial and the effective disorder can grow as the density of the porous 
media decreases. This suggests that instead of changing the pore size a more effective way to control the disorder 
strength is to use mixtures of nematics with other substances whose molecules could "shield" the nematic from the 
random surface. This would decrease the liquid crystal interaction with the porous media. We also showed that at 
certain conditions the porous media can destroy long range order inside a pore. 

We thank Richard Stratt for a useful discussion. This work was supported by the National Science Foundation 
under grant DMR-0131573. 
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